Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid 
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We have used the variational and difTusion quantum Monte Carlo methods to calculate the energy, 
pair correlation function, static structure factor, and momentum density of the ground state of the 
two-dimensional homogeneous electron gas. We have used highly accurate Slater- Jastrow-backflow 
trial wave functions and twist averaging to reduce finite-size effects where applicable. We compare 
our results with others in the literature and construct a local-density-approximation exchange- 
correlation functional for 2D systems. 
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I. INTRODUCTION 

The homogeneous electron gas (HEG) plays a key role in modern condensed-matter theory. It consists of a set of 
electrons moving in a uniform, inert, neutralizing background, and is the simplest fully-interacting quantum many- 
body model of condensed matter. The three-dimensional HEG models the behavior of the conduction electrons in 
metals and semiconductors, while the two-dimensional (2D) HEG models the behavior of electrons confined to surfaces 
and thin layers. Quantum Monte Carlo (QMC) methodsi"^ have long played an important role in establishing the 
ground-state properties of the HEG. Most effort has gone into calculating the ground-state energy of different phases 
as a function of density, in order to establish the zero-temperature phase diagram<^"— In this article, we report QMC 
calculations of some other properties of the 2D HEG of interest to condensed-matter physicists: the pair correlation 
function (PCF), static structure factor (SSF), and momentum density (MD). We also report energy data for high- 
density HEG's. We have confined our attention to the fluid phase, which is the ground state at the densities typically 
encountered in experiments. 

The PCF, especially the contact PCF .g(0), is a key ingredient in generalized-gradient-approximation exchange- 
correlation functionals for density functional theory (DFT) calculations. The PCF has been studied several times 
using QMC^i^i^ii but the value of ^(0) at low densities has proved controversial because electrons approach one 
another infrequently, and the QMC results disagree with values calculated using ladder theoryj^ Our QMC data, 
obtained using a different trial wave function from the earlier calculations, should help to clarify the situation. The 
SSF is related to the PCF by a Fourier transform. SSF data at small k are needed to establish the long-range behavior 
of the PCF. 

The MD of the HEG is of considerable importance in Fermi liquid theory. To our knowledge, the only QMC MD 
data to have been published for the 2D HEG are those of Tanatar and Ceperley;^ which used a relatively simple 
form of trial wave function. (A fit to QMC data generated by Conti^ is reported in Ref. [l^, but no details about the 
calculations are given.) Tanatar and Ceper ley's low-density MD shows a very strange feature: the MD is lower at 
zero momentum than it is at the top of the Fermi edge. It is clearly important to provide new QMC MD data, in 
order to establish whether this is a genuine property of the HEG. 

Finally, we report ground-state energy data for paramagnetic Fermi fluids, which we use to parameterize a local- 
density-approximation exchange-correlation functional for use in DFT studies of 2D systems. 

The rest of this article is arranged as follows. In Sec.|ll]we describe the computational techniques used. In Sec. IIIII 
we present the data we have generated. Finally, we draw our conclusions in Sec. lIVI Densities are given in terms of the 
radius r^ of the circle that contains one electron on average. We use Hartree atomic units {h = |e| = TOe = 47reo = 1) 
throughout this article. All our QMC calculations were performed using the CASINO codeJ^ 

II. QMC CALCULATIONS 
A. Trial wave functions 

In the variational quantum Monte Carlo (VMC) method, expectation values are calculated with respect to an ap- 
proximate trial wave function, the integrals being performed by a Monte Carlo technique. In diffusion quantum Monte 
Carloii^ (DMC) the imaginary-time Schrodinger equation is used to evolve an ensemble of electronic configurations 
towards the ground state. The fermionic symmetry is maintained by the fixed-node approximation)^ in which the 
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Method 


Energy (a.u. / elec.) Var. (a.u. 


) Frac. corr. en. 


HF 


-0.100 222 006 




0% 


SJ-VMC 


-0.148 2110(8) 


0.0196 


96.910(4)% 


SJ-VMC* 


-0.146 80(5) 




94.1(1)% 


SJB-VMC 


-0.149 3851(6) 


0.007 74 


99.282(4)% 


SJB-VMC* 


-0.148 80(5) 




98 1(1)% 


SJ-DMC 


-0.149177(8) 




98.86(2)% 


SJ-DMC* 


-0.149 00(5) 




98.5(1)% 


SJ-DMCt 


-0.149 134(9) 




98.77(2)% 


SJB-DMC 


-0.149 741(2) 




100% 


SJB-DMC* 


-0.149 55(5) 




99.6(1)% 


SJB-DMCt 


-0.149 518(9) 




99.55(2)% 



TABLE I: Energy, variance, and percentage of correlation energy retrieved using different methods for a 58-electron param- 
agnetic Fermi fluid of density parameter = 5 a.u. Twist averaging has not been used. "HF," "SJ-VMC," "SJB-VMC," 
"SJ-DMC," and "SJB-DMC" stand for Hartree-Fock theory, VMC with a Slater- Jastrow wave function, VMC with a Slater- 
Jastrow-backflow wave function, DMC with a Slater-Jastrow wave function, and DMC with a Slater- Jastrow-backflow wave 
function, respectively. The DMC energy data have been extrapolated to zero time step. The data marked with an asterisk 
were produced by Kwon et al.^ while the data marked with a dagger were generated by Attaccalite et al^ at a time step of 
0.1 a.u. (i.e. their data were not extrapolated to zero time step). The fraction of the correlation energy retrieved is computed 
on the assumption that our Slater-Jastrow-backflow DMC calculation retrieves 100% of the correlation energy. 



nodal surface of the wave function is constrained to equal that of a trial wave function. The VMC algorithm gener- 
ates electron configurations distributed according to the square of the trial wave function, while the DMC algorithm 
generates configurations distributed as the product of the trial wave function and its ground-state component. 

Our trial wave functions consisted of Slater determinants of plane- wave orbitals multiplied by a Jastrow correlation 
factor. The Jastrow factor contained polynomial and plane-wave expansions in electron-electron separations^ The 
orbitals in the Slater wave function were evaluated at quasiparticlc coordinates related to the actual electron positions 
by backflow functions consisting of polynomial expansions in electron-electron separations^ The wave functions were 
optimized by variance minimizatiorS^S^ and linear-least-squares energy minimizations^ 

We simulated HEG's in finite, square cells subject to periodic boundary conditions. The many-body Bloch theorenS^ 
states that the wave function satisfies 

^(ri, . . . , r, + R,, . . . , r^) cxp(ik, • Rs)4'(ri, . . . , tat), (1) 

where Rs is a simulation-cell lattice vector and kg is the simulation-cell Bloch vector. In some of our calculations, 
and in previous QMC studies of the 2D HEGj^i^S^ it has been assumed that = 0. However, in our calculations 
of the energy, PCF, and SSF we performed twist averaging, in which expectation values are averaged over k^ in the 
first Brillouin zone of the simulation celli^ This procedure greatly reduces single-particle finite-size errors caused by 
shell-filling effects. 

The high quality of our trial wave functions is demonstrated in Table HI which shows QMC energies achieved using 
different levels of wave function for a 58-electron paramagnetic Fermi fluid of density parameter = 5 a.u. Backflow 
functions change the nodal surface of the trial wave function and can therefore improve the fixed-node DMC energy. 
In practice we find that backflow lowers the DMC energy substantially. Our VMC energies are significantly lower than 
those of Kwon et al.r^ as is our Slater-Jastrow-backflow DMC energy. On the other hand, the Slater-Jastrow-backflow 
DMC energy of Attaccalite et al^ is higher than that of Kwon et al. Our Slater-Jastrow DMC energy is slightly lower 
than that of Attaccalite et ai, which in turn is lower than that of Kwon et al. Since the nodal surface is the same in 
the three calculations, these DMC energies really ought to agree. However, the trial wave function used by Kwon et 
al. is very much poorer than ours, as can be seen by comparing the VMC energies. Time-step and population-control 
biases in their DMC energies must be much greater, which may explain the discrepancy. The results of Attaccalite 
et al. have not been extrapolated to zero time step, hindering comparison. The VMC and DMC results of Rapisarda 
and Senatore^ are in very close agreement with those of Kwon et alJ^ 

We have optimized a three-electron term in the Jastrow factor (together with the two-electron Jastrow terms and 
backflow functions) for a paramagnetic 58-electron HEG at rs = 5 a.u. The three-electron term lowered the non- 
twist- averaged VMC energy from —0.1493851(6) to —0.1495111(5) a.u. per electron. The DMC energies at a time 
step of 0.1 a.u. without and with the three-body Jastrow factor are —0.149742(2) and —0.149740(2) a.u. per electron. 
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respectively. As expected, the inclusion of the three-body term makes an insignificant difference to the DMC energy, 
because the DMC energy depends only on the nodal surface of the trial wave function, which is not directly affected 
by the Jastrow factor. We have therefore not used three-electron terms in our production calculations. 



B. Evaluating expectation values 

1. Evaluating the MD 

Let ^'(R) be the trial many-electron wave function, where R = (ri, . . . ,rjv). Suppose the first electrons are 
spin-up and the remainder are spin-down. The MD of spin-up electrons can be evaluated as 

where the angled brackets denote an average over the set of electron configurations generated in the VMC and DMC 
algorithms (which arc distributed as and l^'fE'oli respectively, where $o is the ground-state component of '5). We 
have restricted our attention to paramagnetic and fully ferromagnetic HEG's, so the total MD is equal to the spin-up 
MD. The integral in the expectation value of Eq. ([2]) is estimated by Monte Carlo sampling at each configuration R 
generated by the QMC algorithms, and the results are averaged. The use of a finite number of points in the evaluation 
of the integral at each R does not bias the QMC estimate of p(k). 

Suppose our finite simulation cell has area A and the simulation-cell Bloch vector is kj,. We may write 

«'(r,r2, . . . ,rAr) 1^ ,i ^ i /Q^ 
= -2^CG(R)exp[z(G + k,) -r], (3) 

where the {G} are the simulation-cell reciprocal lattice points. Hence it is clear that p{\i) is only nonzero if k = G + kg 
for some G. The MD is only defined for a discrete set of momenta at any given kg. One cannot twist average as such; 
instead, altering ks leads to the MD being defined at a different set of momenta. We simply report MD's obtained 
using ks = (i.e., no twist was applied). 



2. Evaluating the SSF 



The SSF may be evaluated as 

Sik) = 1 [(n(k)n(-k)) - (n(k)) (n(-k))] , (4) 

where 

n(k) = ^exp(-«k • ri) (5) 

i 

is the Fourier transform of the density operator. S'(k) is only nonzero at simulation-cell G vectors, even if the 
simulation-cell Bloch vector is nonzero. We can twist average when we calculate 5(G). 



3. Evaluating the PCF 

The spherically averaged PCF is 



5W = i^(E'^(l-'-r.|-'^)), (6) 



which can be evaluated by binning the electron-electron distances in the configurations generated by the QMC 
algorithms. Twist averaging introduces no complications. 
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FIG. 1: (Color online) Twist-averaged DMC energy against time step for paramagnetic Fermi fluids of density parameter 
rs = 5 a.u. at different system sizes A^. "SJ" and "SJB" refer to Slater-Jastrow and Slater- Jastrow-backflow wave functions, 
respectively. 



4- Extrapolated estimation 

If A is an operator that docs not commute with the Hamiltonian then the errors in the VMC and DMC estimates 
^VMC s-nd Ay)mc oi the expectation value of A are linear in the error in the trial wave function; however, the error 
in the extrapolated estimate 2Admc ^ ^vmc is quadratic in the error in the trial wave functioni^ We have used 
extrapolated estimation in most of our calculations of expectation values. Examples of extrapolation are shown in 
Figs. [5]and[8l and the upper panel of Fig. S) In each case the VMC, DMC, and extrapolated estimates are in good 
agreement, implying that the error resulting from the use of a DMC mixed estimate is small. Gori-Giorgi et ali^ used 
reptation2^ QMC to accumulate the PCF and SSF, in which pure expectation values are obtained with respect to the 
fixed-node ground-state wave function, so that extrapolation is unnecessary. 



C. Time-step and population-control biases 



Finite-time-step errors in the twist-averaged DMC energy were removed by linear extrapolation to zero time step. 
An example is shown in Fig. [T] it can be seen that the time-step bias is in fact very small in any case. We checked 
that the other expectation values were converged with respect to the time step: see Figs. |3]and[ni We used a target 
population of 1600 configurations in all our DMC calculations, making population-control bias negligible. 



D. Finite-size bias 



Expectation values obtained in a finite A'^-electron cell subject to periodic boundary conditions differ from the 
corresponding infinite-system values because of "single-particle" shell-filling effects, as well as the neglect of long- 
ranged correlations and the compression of the exchange-correlation hole into the simulation cell. Single-particle 
finite-size effects can be removed from the energy, the SSF, and the PCF by twist averaging, as explained in Sec. Ill Al 
We have recently demonstrated that the finite-size error in the energy per particle in a 2D HEG falls off as N~^/^, 
enabling us accurately to extrapolate QMC energies to infinite system size<^ For the PCF, SSF, and MD we simply 
verified that the QMC data had converged with respect to system size (see Figs. [S] and [51 and the lower panel of Fig. 

a. 
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rs (a.u.) A*' DMC energy (a.u. / elec.) 



1 


50 


-0.212 5(2) 


1 


74 


—0.212 2(2) 


1 


114 


-0.211 6(3) 


1 


oo 


-0.210 4(6) 


5 


58 


—0.149 95(2) 


5 


114 


—0 149 77(1) 


5 


oo 


-0.149 63(3) 


10 


58 


-0.085 504(5) 


10 


74 


-0.085 52(1) 


10 


114 


-0.085 445(3) 


10 


oo 


-0.085 399(6) 



TABLE II: Twist-averaged DMC energy, extrapolated to zero time step, for A''-electron paramagnetic Fermi fluids of density 
parameter r^. Where — oo, the DMC energy has been extrapolated to infinite system size. 



III. RESULTS 
A. Energies 

DMC energies of paramagnetic Fermi fluids at different densities and system sizes are shown in Table HT] Our 
results for the energies of different phases of the 2D HEG at low density are reported elsewhere.— At = 5 and 10 
a.u., Rapisarda and Senators^ obtained infinite-system energies of —0.1490(1) and —0.08512(2) a.u. per particle using 
DMC with a Slater-Jastrow wave function. Kwon et alt^ obtained DMC energies of —0.2098(3), —0.1495(1), and 
—0.08536(2) a.u. per electron at rg = 1, 5, and 10, respectively, using a Slater- Jastrow-backflow wave function. Our 
DMC energies are somewhat lower than these data, as expected from the results shown in Table HI 

Let the correlation energy per electron Ec be the difference between the ground-state energy per electron and the 
Hartree-Fock energy. We fit the form proposed by Attaccalite et al^ to our correlation energies for paramagnetic 
EEC's: 

E,^Ao + (Bor, + Corl + Dorl) log ( 1 + — ^ ) , (7) 

where Aq = —0.1925, Bq = ■\/2(10 — 37r)/(37r), and Dq = —AqHq. We fit to the infinite-system DMC energies shown 
in Table inland also to the DMC energies of low-density paramagnetic HEG's reported in ReL [B] (at rg = 20, 25, 30, 
35, and 40 a.u.). Our fitting parameters are shown in Table Hill and the correlation energies of paramagnetic Fermi 
fiuids obtained by different authors relative to that of Attaccalite et al. are shown in Fig. 12] Our correlation energies 
are lower than those of the other authors because of our use of flexible backflow functions. Equation (O can be used 
as a local-density-approximation exchange-correlation functional in DFT calculations for 2D systems. 

Unlike Attaccalite et al.r^ we fit Eq. (|71) to paramagnetic data only; we do not attempt to calculate the spin- 
polarization-dependence of the energy of the HEG. 



B. MD's 



The MD's of paramagnetic HEG's are shown in Fig. 121 and a more detailed graph of the MD at = 5 a.u. is shown 
in Fig. m The upper panel of Fig. SI demonstrates that the extrapolated estimate is accurate and that the DMC 
results are converged with respect to the time step. It is clear from the lower panel of Fig. 21 that, although backflow 
makes a significant improvement to the QMC energy estimates, it has very little effect on the MD. The inclusion of 
backflow results in a small transfer of weight to wavevectors above the Fermi edge, as expected, because a greater 
fraction of correlation energy is retrieved. It can also be seen that the MD's obtained at different system sizes are in 
agreement. We have therefore plotted data obtained at different system sizes together in Fig. I3l To our knowledge, 
the only previous QMC studies of the MD of the 2D HEG are those of Tanatar and Ceperley^ and Conti<^ At = 10 
a.u., Tanatar and Ceperley found that the MD at small wave vectors is lower than the value near the Fermi edge. 
Our data do not show this unusual feature. Tanatar and Ceperley used a relatively inflexible Slater-Jastrow wave 
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Parameter 


Value 


^0 


-0.192 5 


Bo 


0.086 313 631 


Co 


0.069 795 68 


Do 


0.0 


Eo 


1.053 100 3 


Fo 


0.040 69122 


Go 


0.360 595 3 


Ho 


0.0 



TABLE III: Fitting parameters in Eq. (O for the correlation energy of a paramagnetic HEG. Equation was fitted to the 
infinite-system DMC energies given in Table [III and the low-density energy data given in Ref.[^. It was found that fixing Hq = 
did not affect the quality of the fit. 



Attaccalite et al. 

• Present work 

♦ Kwon et al. 

■ Rapisarda & Senatore 
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FIG. 2: (Color online) Correlation energy of a paramagnetic Fermi fluid relative to the results of Attaccalite et al^ The results 
obtained by Kwon et ali^ and Rapisarda and Senators^ are shown for comparison. The results of Tanatar and Ceperlejj^ are 
not shown, because they are systematically too low in energy.— 



function, which may be the reason for the discrepancy. Giuliani and Vignale^^ quote a formula for the MD, which 
was obtained by fittingto QMC data generated by Contii^ We have fitted our MD's to a simplified version of the 
form suggested in Ref. Iiol : 



4g(0)r~ 



a2X~ 



+ (ay + a%X + OgX' 



040; 
') exp 



ii x<V2 



(8) 



where x = Vgk and the Ui are fitting parameters. g{Q) is the contact PCF, which we evaluated using Eq. (|9]). The 
fitted parameters are give n in Table IIVI Our values for the discontinuity at the Fermi edge are slightly smaller than 
those reported in Ref. [Tfl- 



C. SSF's 



VMC and DMC SSF's of a 50-electron paramagnetic HEG at = 1 a.u. are shown in Fig. [5] It can be seen that 
the difference between the VMC and DMC data is in most cases smaller than the difference between the two sets 
of DMC data, implying that the errors due to extrapolated estimation are small. On the other hand, the difference 
between the data with kg = and the twist-averaged data is significant. In particular, the former has some unusual 
features close to integer multiples of the Fermi wave vector, one of which is shown in the inset to Fig. [5l Elsewhere, 
twist averaging has only a small effect. At all densities we find that the twist-averaged SSF in a 50- or 58-electron 
cell is in agreement with the SSF in a 114-electron cell, as can be seen in Fig. \5\ Since the statistical errors are less 
significant at the smaller system sizes, we have used A'^ = 50 or 58 electrons in the high-density data reported below. 
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FIG. 3: (Color online) Extrapolated MD p{k) for paramagnetic Fermi fluids, fcp = V^/ts is the Fermi wave vector of the 
paramagnetic fluid and pp ~ r^/ (27r) is the value of the Fermi distribution. The results were obtained using a Slater- Jastrow- 
backflow wave function and a variety of system sizes with A'' > 50 in each case. Twist averaging was not used. For comparison, 
we have plotted the MD obtained by Tanatar and Ceperlej«^ (open symbols) and Eq. (8.133) of Ref. [l3 (dotted lines). 
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FIG. 4: (Color online) MD p{k) for paramagnetic Fermi fluids at Ts = 5 a.u., obtained using VMC and DMC at different system 
sizes A''. Backflow was used, except for the data labeled "SJ." The upper panel shows the results obtained at A'^ = 58 using 
VMC and DMC with different time steps r. The lower panel shows the effect of changing the system size and using backflow 
on the extrapolated MD. Twist averaging was not used, fcp = V2/rs is the Fermi wave vector of the paramagnetic fluid and 
pp — r1/{2-K) is the value of the Fermi distribution. 
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rs (a.u.) 


1 


5 


10 


30 


ao 


1.950 


1.649 


1.410 


0.974 5 


ai 


-0.073 42 


-0.038 99 


0.336 6 


-0.053 5 


a2 


0.280 5 


0.07418 


-1.199 


0.150 9 


as 


-0.388 4 


-0.192 


1.148 


-0.388 8 


a4 


0.136 5 


0.02198 


-0.436 3 


0.147 3 


as 


1.171 


1.017 


1.035 


1.379 


a? 


0.164 8 


1.682 


2.091 


1.428 


ag 


-0.113 5 


-1.282 


-1.566 


-0.826 4 


ag 


0.021 19 


0.277 3 


0.343 8 


0.1599 


Z 


0.866 


0.398 


0.209 


0.055 5 



TABLE IV: Fitting parameters in Eq. ^ for the MD and discontinuity Z in the MD at the Fermi edge. The MD is normahzed 
such that the Fermi distribution is 1. 
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FIG. 5: (Color onhne) SSF S{k) for a paramagnetic HEG of density parameter 7-3 = 1 a.u. fcp = is the Fermi wave 

vector of the paramagnetic fluid. The system size was A'' = 50 electrons, except where indicated otherwise. Results obtained 
with ks = ( "PBC" ) and twist averaging ( "TA" ) are shown. A Slater- Jastrow-backflow wave function was used in each case. 
The inset shows in greater detail one of the regions in which the twist-averaged and non-twist-averaged data disagree. 



The SSF's of paramagnetic and ferromagnetic fluids at low density are shown in Fig. [6l while the SSF's of para- 
magnetic fluids at high density are shown in Fig. [T] It can be seen that a peak in the SSF at about k = 2.62fcF 
appears at low density, perhaps due to incipient Wigner crystallization. Our SSF's are in good agreement with those 
of Gori-Giorgi et alX 



D. PCF's 



We compare PCF's obtained at different system sizes using different QMC methods for a paramagnetic HEG at 

= 5 a.u. in Fig. [5] The difference between the VMC and DMC PCF's, and the difference between extrapolated 
PCF's obtained with and without backflow is small, implying that the error in the extrapolated PCF is small. Twist 
averaging has a small effect on the PCF, but the twist-averaged PCF's at = 58 and N = 114 electrons are very 
similar, implying that the finite-size error in the twist-averaged PCF at = 58 is small. We have also verified that 
the PCF is converged with respect to the DMC time step. 

Our PCF's are shown in Figs. [S] (high density) and [TUl (low density), along with the results of Gori-Giorgi et alJ- 
Our PCF's are in good agreement with those of Gori-Giorgi et al. (as expected, from the SSF results in Fig. [7]). 

The contact PCF g(0) is especially important in the construction of gcneralized-gradient-approximation exchange- 
correlation functionalsji^ We give our g{0) values in Table |Vl and we plot rsg{0) against in Fig. [TTJ along with 
some other results in the literature. Because our PCF's have converged with respect to system size, and our VMC and 
DMC results agree with each other when backflow is used, we have simply averaged our VMC and DMC g{0) data at 
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FIG. 6; (Color online) Extrapolated SSF S{k) for HEG's of density parameter Ts and spin polarization C,. fcp = \/2/rs is the 
Fermi wave vector of the paramagnetic fluid. The results were obtained using a Slater- Jastrow-backflow wave function and 
twist averaging. System sizes of A*' = 90, 114, 90, and 114 were used in the paramagnetic calculations at = 20, 25, 30, and 
35 a.u., respectively, and a system size of A'^ = 45 was used in the ferromagnetic calculation at rs = 20 a.u. The inset shows 
the peak in greater detail. 
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FIG. 7: (Color online) Extrapolated SSF S{k) for paramagnetic HEG's at high density, fcp = \/2/rs is the Fermi wave vector. 
Slater- Jastrow-backflow wave functions and twist averaging were used. The system sizes are = 50, 58, and 58 at rs = 1, 5, 
and 10 a.u., respectively. The curve marked "HF" shows the Hartree-Fock SSF. The dotted lines show the SSF obtained by 
taking the Fourier transform of the PCF data of Gori-Giorgi et ali^ 



different system sizes in order to reduce the statistical noise in our estimate of the contact PCF. Our results are in 
reasonably good agreement with the fit to the earlier QMC data of Gori-Giorgi et alJ- and also with the expression for 
g{0) obtained using ladder theory by Nagano et al^ Interestingly, our results clearly disagree with the more recent 
calculation of g(0) within ladder theory by Qian,— which involved fewer approximations than the work of Nagano et 
al. The close agreement between our results and those of Nagano et al. must therefore be regarded as a coincidence. 
The fact that our QMC calculations, using a different trial wave function, are consistent with the data of Gori-Giorgi 
et al. strongly suggests that the QMC results for ^(0) are reliable, whereas ladder theory is of limited use at low 
densities. Our results are also in clear disagreement with the formula proposed by Polini et al.^ which interpolates 
between the results of ladder theory at high-density (where it should be exact) and a partial-wave analysis at low 
density. 

The fit to our g{{)) data shown in Fig. [TT] is 

2. \l + ars + brj.+ cr^l <l 
where A = -0.25724, B = 0.071116, and E = 0.98553 were obtained by fitting. Polini et al^ have shown that 
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FIG. 8: (Color online) VMC, DMC, and extrapolated ("Ext.") PCF's g{r) for a paramagnetic Fermi fluid of density parameter 
rs = 5 a.u. at different system sizes A*'. Twist averaging was used in the curves labeled "TA," but not in the one labeled "PBC." 
Slater-Jastrow and Slater-Jastrow-backfiow wave functions were used in the curves labeled "SJ" and "SJB," respectively, 
fcp = V^/tb is the Fermi wave vector of the paramagnetic fluid. 
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FIG. 9: (Color online) Extrapolated total PCF g{r) for paramagnetic Fermi fluids of density parameter r^. kp — yf^jrs is the 
Fermi wave vector of the paramagnetic fluid. "HF" stands for Hartree-Fock theory. Twist averaging was used and the QMC 
calculations were performed at system sizes of A'^ = 50, 58, and 58 electrons at = 1, 5, and 10 a.u. The dotted lines show 
the data of Gori-Giorgi et al.j^ which are almost indistinguishable from our data. 



linir^_^o .9(0) = (1/2) [1 - 1.372rs], so we have set a = -1.372 and determined b = 0.997618888 and c = -0.3218467056 
by matching the value and derivative of 17(0) at = 1 a.ui^ 

IV. CONCLUSIONS 

We have studied the ground-state properties of the fluid phases of the 2D HEG using QMC. We used highly 
accurate trial wave functions and dealt with finite-size effects by twist averaging. Twist averaging removes some 
strange features in the SSF, but our PCF's and SSF's are in good agreement with analytic fits to earlier QMC 
data^ confirming the accuracy of these formulas. Our MD's show some qualitative differences from earlier QMC 
results;^ however; in particular, we do not observe an increase in the MD as the Fermi edge is approached at low 
density. Finally, we have reported DMC energy data for the high-density 2D HEG, which wc used to construct a new 
exchange-correlation functional for 2D DFT calculations. 
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FIG. 10: (Color online) Extrapolated total PCF g{r) for Fermi fluids of density parameter Vs and spin polarization C,. fcp — 
V2/ri, is the Fermi wave vector of the paramagnetic fluid. "GMB" denotes the work of Gori-Giorgi et ali- (shown by dotted 
lines), while "HF" stands for Hartree-Fock theory. Twist averaging was used and the QMC calculations were performed at 
system sizes of A'' = 90, 114, 90, and 114 at = 20, 25, 30, and 35, respectively. 



(a.u.) (?(0) 



1 
3 
5 
7 
10 



0.1517(4) 
0.022 7(2) 
0.005 0(3) 
0.00125(5) 
0.000 24(2) 



TABLE V: Contact PCF of paramagnetic HEG's at five different densities. Twist-averaged VMC and DMC results obtained 
at different system sizes were averaged to obtain these data. 
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FIG. 11: (Color online) Contact PCF glfi) against density parameter r^, as calculated by different authors: the present work 
(see Table |V}, Gori-Giorgi et al.^ Nagano et al.^ Polini et al.^ and Qian.— 
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